Receptor model-based sources and risks appraisal of potentially toxic elements in the urban soils of Bangladesh

Rapid urbanization and industrial development have prompted potentially toxic elements (PTEs) in urban soil in Bangladesh, which is a great concern for ecological and public health matters. The present study explored the receptor-based sources, probable human health and ecological risks of PTEs (As, Cd, Pb, Cr, Ni, and Cu) in the urban soil of the Jashore district, Bangladesh. The USEPA modified method 3050B and atomic absorption spectrophotometers were used to digest and evaluate the PTEs concentration in 71 soil samples collected from eleven different land use areas, respectively. The concentration ranges of As, Cd, Pb, Cr, Ni, and Cu in the studied soils were 1.8—18.09, 0.1—3.58, 0.4—113.26, 0.9—72.09, 2.1—68.23, and 3.82—212.57 mg/kg, respectively. The contamination factor (CF), pollution load index (PLI), and enrichment factor (EF) were applied to evaluate the ecological risk posed by PTEs in soils. Soil quality evaluation indices showed that Cd was a great contributor to soil pollution. The PLI values range was 0.48–2.82, indicating base levels to continuous soil degradation. The positive matrix factorization (PMF) model showed that As (50.3 %), Cd (38.8 %), Cu (64.7 %), Pb (81.8 %) and Ni (47.2 %) were derived from industrial sources and mixed anthropogenic sources, while Cr (78.1 %) from natural sources. The highest contamination was found in the metal workshop, followed by the industrial area, and brick filed site. Soil from all land use types revealed moderate to high ecological risk after evaluating probable ecological risks, and the descending order of single metal potential ecological risk was Cd > As > Pb > Cu > Ni > Cr. Ingestion was the primary route of exposure to potentially toxic elements for both adults and children from the study area soil. The overall non-cancer risk to human health is caused by PTEs for children (HI=0.65 ± 0.1) and adults (HI=0.09 ± 0.03) under USEPA safe limit (HI>1), while the cancer risks from exclusively ingesting As through soil were 2.10E-03 and 2.74E-04 for children and adults, respectively, exceeding the USEPA acceptable standard (>1E-04).


Introduction
Soil pollution by PTEs is a matter of great concern because of their source abundance, richness, slow breakdown, higher accumulating rate, and probable eco-toxic nature for humans and other biota [1][2][3]. Soil is a valuable natural reserve for the existence of biotic components [2,4], which is being contaminated by PTEs from manmade activities, including waste disposal, sewage effluent, agricultural practices, industrial activities, vehicular exhaust, etc [5]. In comparison to water and air, soil can acquire larger amounts of organic and inorganic substances because of their diverse influencing factors and increasing area [6]. Nowadays, soil pollution by PTEs in urban areas is a highly remarkable environmental issue for developed and developing countries [7]. In urban areas, soil pollution is gradually increasing because of rapid urbanization and industrialization activities [8]. This problem is common in Bangladesh because many landowners switch their agricultural land for diverse purposes, such as shopping malls, markets, workshops, storage sites, etc., for high economic benefits [9]. The distribution of toxic metals in urban area soil varies from place to place due to the characteristics of parent materials, climate, soil morphology, and anthropogenic interferences [10,11]. Therefore, urban soil act as a sink of pollution due to the continuous receiving of an excessive amount of PTEs, consequently leading to the changing in soil chemistry, soil ecology, and soil functioning capacity [12][13][14]. The studies of probable ecological risk evaluation and Source attribution of PTEs in contaminated soil have attained more attention [15]. Diverse environmental risk evaluation approaches apply when data is taken from non-polluted areas [2,8,[16][17][18]. Soil serves as a range of habitats, stores massive contaminants, and significantly contributes to environmental and geochemical speciation, so surface soil sampling is a fast, inexpensive, effective, and significant information-carrying system for evaluating the quality of terrestrial ecosystems [19,20]. Conversely, the bioavailability of toxic metals in urban soil can pose substantial human health risks if they entered the human body via dermal contact, incidental ingestion, and inhalation [21][22][23]. The PTEs exposure level and potential health risk depend on different toxicity influencing factors, such as chemical type, age, gender, dose-time relation, mode of exposure, metabolic rate and distribution within the body of the host etc. [24]. Children and infants are more vulnerable to PTEs exposure and potential health risk than adults because they may inhale, ingest and dermally absorb greater amounts of heavy metals by involving their hands, mouth, or even the whole body during their outdoor and play activities in the urban environment [25][26][27]. After exposure, PTEs are stored in the human body (eg. bloodstream, gastrointestinal and intestinal membranes) that either create heavy metal-associated health hazards or act as a co-factor of several diseases, while older people are highly threatened due to their lower disease defence capacity [3,28,29]. So, to determine and define the qualitative (cancer-causing and non-cancer-causing health concerns) and quantitative (mapping and calculation) risk levels, methods for risk assessment approach are required [17,26]. Moreover, the long-period accumulation of PTEs not only degrades the soil quality and their adjacent environment but also poses a significant adverse effect on human health. Therefore, it is essential to investigate the distribution of PTEs and the identification of sources in different types of urban soil usage to evaluate the level of contamination and develop remediation strategies. Conversely, detailed source dispersal of PTEs in urban land use soil is far away from scientific consideration yet this study is vital for public health concerns. Due to the variety and complexity of the pollution sources, using a single scientific method to draw pollution zones is challenging. To identify the specific and acceptable results of PTEs sources, this study applied principle component analysis (PCA), and receptor model (e.g. PMF). This study's endeavour to incorporate different receptor models (PCA, and PMF) may be an appropriate strategy for the different pollution sources scattered throughout Bangladesh urban soils. The Jashore district was selected for this study as the location that has drawn the most attention owing to experiencing a variety of environmental issues as a result of the fast expansion of infrastructure and industrial activity. Monitoring the distribution of poisonous metal levels in urban soil surfaces is therefore required to compare with non-polluted background values and evaluate the environmental health hazard and ecological sustainability for this area.
Though several studies have been conducted in agricultural soils of Bangladesh for determining the metal concentration such as industrial areas in Dhaka [30], Dhaka Export Processing Zone [31], Dhaka Export Processing Zone Area [32], Patuakhali district [5], industries area of Dhaka city [33], the evidence on PTEs in the urban soil by land-use changes in Jashore district especially using PMF model is inadequate. Therefore, the objectives of this study aimed to identify (1) the concentrations of PTEs (As, Cd, Pb, Cr, Ni, and Cu) in different land-use urban soils; (2) to find out the potential sources of PTEs in the urban soil using receptor model-based sources; and (3) to analyze the possible risks of PTEs to ecological and human health.

The study area
Jashore is one of the more significant districts of the Khulna division, located in Bangladesh southwestern region. Geographically it is located at 22 • 10 ´and 22 • 28 ´east longitudes and between 89 • 16 ´and 89 • 64 ń orth latitudes (Fig. 1). This study area was exposed to a variety of pollutants from the outflow of urban areas, heavy traffic, and airborne fallout, as well as solid waste and untreated sewage from numerous industrial sites (such as brickfields, chemical industries, tanneries, and, textile). The country's biggest Land port (Benapol) is situated in this area, mainly used for import and export trading between Bangladesh and India. Lithologically, this area is formed by the Ganges deltaic action, and it is also situated at Gopalganj-Khulna Bills agro-ecological zone no 14.
Between the Ganges River Floodplain as well as the Ganges Tidal Floodplain, the region comprises large low-lying regions, and clay loam soils are predominant in the Jashore district [34]. The yearly rainfall is 1537 millimetres (60.5 in.) and the yearly mean temperature ranges from 15.4 • to 34.6 • C (59.7-94.3 • F) [35], which significantly influences the concentration and dispersal of metal and metalloids in elution from both soil and surface runoffs [36]. Conversely, the whole area of the district is 2606.98 square km, 7 m elevated from sea level, with a population of a total of 2764,547 [37], and the main rivers are the Bhairab, Chitra, Betna, Kopotakkho, and the Mukteshwari, while this district consists of eight Upazila [38].

Collection, processing, and metal extraction of soil samples
This research was performed in 11 different land uses sites i.e., brickfield (BF), playground (PG), market (M), park area (PA), industrial area (IA), petrol pump (PP), metal workshop (MW), bus stand (BS), power station (PS), residential area (RA), and waste dumping area (WDA), where surface soil samples were collected during the January to March 2021 from 71 sampling points (Fig. 1), at a depth of 0-15 cm using a stainless-steel scoop. Three subsoil samples were taken at each location following the standard protocol [39]. These sub-soil samples were appropriately mixed to make a composite sample. The study soil sample was taken from these composites and stored in nonmetallic sample collecting bags to avoid unwanted contaminations. The obtained samples were dried for 24 h at 80 • C in the oven (Labtech LDO-150 F, Korea) then crushed using a crusher and sieved (2 mm) [40]. For lab analysis, the samples were then placed in Pyrex-sealed glass vials. The pre-industrial soil sample for background value was collected and prepared by following the procedures of Islam et al., 2015. Analytical grade chemicals (Merck, Germany) were used throughout the whole analysis. Modified USEPA Method 3050B [41] was used to carry out the entire soil digestion to determine the presence of toxic metals. To summarize, 2.0 g of every soil sample was sequentially digested utilizing concentrated 10 mL HNO 3 (69 %, Merck, Germany; 95 0 C for 2 h), 5 mL H 2 SO 4 (98 %, Merck, Germany; 250 0 C for 2 h), 5 mL HClO 4 (70 %, Merck, Germany; 300 0 C for 2 h) and 3 mL H 2 O 2 (30 %, Merck, Germany; 300 0 C for 30 min). Following digestion, solutions were extracted using Whatman no. 41 filter paper (pre-washed with 0.1 M HNO 3 ), and the final volume was created using 100 mL of double-distilled water.

Chemical investigation and quality control
The physicochemical composition of soil samples' such as electrical conductivity (EC) (HACH Sension − 156; multi-parameter, USA), soil texture (sand, silt, and clay), pH meter (Model: MARTINI instruments, pH 56 pHWP, USA), percentage of organic carbon (OC), and percentage of organic matter (OM) was determined using standard procedure, adopted from Kormoker et al. [17]. A hydride Generator, graphite furnace and air-acetylene flame incorporated

Appraising the contamination level of PTEs
Different pollution evaluation indices, including the enrichment factor (EF), geo-accumulation index (I geo ), contamination factor (CF i ), degree of contamination (C d ), pollution load index (PLI), and, potential ecological risk evaluation index (PERI) were used in this study to assess the degree of pollution, details presented in Table 2 and pollution/risk evaluation category are given in Table S2.

Human health risk assessment
This study estimates the human health risks (non-carcinogenic and carcinogenic) due to exposure to contaminated soil via ingestion, dermal contact, and inhalation exposure pathways for children and adults. USEPA [39,48] proposed health risk assessment method was applied ( Table 2), and model parameters and references dose are presented in Table 3 and Table S3, respectively.

Positive matrix factorization (PMF)
PMF (EPA PMF version 5.0) is a mathematical receptor model, used in this study to distribute the source of toxic metals in the urban soil. It implements well to calculate summary and source contribution according to reliable factorization algorithms; mathematically it can be stated as Eq. (18): (18) Where, x iy is the i th species value calculated in the y th sample; f ay is the input of a th source to y th sample; g ia is the value of i th species from y th source, and e iy is the methodical error. The objective of this model was to discover the values for g ia and f ay , which best repeat the measured value x iy . These data were perfected until the lowest Q value was attained, where Q is defined as Eq. (19).
Here the α iy states the "uncertainty" in the i th species of sample number y. heavy metals concentration and their agreeing uncertainties are the input data for the model, which are evaluated based on PMF fundamentals (version 5.0) and User Guide. The data below the detection limit (MDL) are replaced with the methods as follows in Eqs. (20)(21).
Once the lower Q values, higher R 2 values, and clear interpretability of the factors were attained, the data were run 20 times using this model with a random seed and a varied number of factors (varying from 3 to 6). All runs' Q (robust) and Q (True) values were well examined, and run 13 was chosen for factor extraction, with steps beginning at least 127. Positive matrix factorization model accuracy data are presented in Table S4.

Statistical analysis
Multivariate statistical techniques such as principal component analysis (PCA), cluster analysis (CA), and Pearson's correlation matrix (PCM) were used to determine the most likely sources of toxic harmful metals in the soil. To explain the study findings, PCA was performed Table 1 Pollution evaluation indices.

Evaluation indices Equation Explanation
Reference (C x / M Fe ) sample = the ratio of the concentration of toxic metals (C x ) to that of iron (M Fe ) in the soil sample (C M /C Fe ) background = the same reference ratio in the pre-industrial sample Chakraborty et al. [42], Islam et al. [7], Sutherland [43] Geo-accumulation index (I geo ) (2) C xi = the studied toxic metal concentration of the i th parameter B xi = the geochemical background value of the toxic metal of the i th parameter in the preindustrial soil samples of the study area Muller [44], Islam et al. [45] Contamination factor Hakanson [46] Degree of contamination (C d ) Hakanson [46] The pollution load index (PLI) n = number of studied metals Tomlinson et al. [47] Potential ecological risk evaluation index (PERI)

Physiochemical characteristics and concentration of PTEs in soil samples
The soil pH values varied from 7.15 to 8.36 (Table S5), indicating that the slightly alkaline might be the decaying of organic matter as well as discharging of wastes, effluents, chemicals and salt etc in soil [51]. Soil Resources Development Institute (SRDI), Bangladesh, suggested that this study area's soils are in the non-saline to low saline category. The percentage of organic matter ranged from 1.09 to 5.65% (Table S3), due to the deposition or decomposition of the larger amount of waste and sewage sludge. Several studies recommended that polluted soil contains more organic matter than non-polluted soil [51]. Based on the USDA soil texture classification, this study area's soil is sandy loam ( Table S5). The PTEs concentration ranges in the soil samples of different land uses of the study areas were as follows: As  Table 2 Human health risks assessment methods via ingestion, dermal contact and inhalation exposure routes.
(3.82-212.57 mg/kg) ( Table 1.). Changes in the point and non-point sources of pollution may be responsible for the fluctuation of PTEs. The descending order of the mean PTEs concentration in soil samples was Cu > Cr > Pb > Ni > As > Cd (Table 4). The relatively high standard deviation for some metals (eg. Cr, Cu, and Pb) in different land uses exhibited relatively substantial spatial variability in these metals [52]. Potential toxic element concentrations in study locations followed the decreasing order of site MW Various metallic waste, burning of fossil fuels, emissions, and waste from brickfields and metal workshop areas are the main sources of As [55].
Higher concentrations of Cd found in soil was 2.05 ± 0.84, 1.25 ± 0.86, 1.24 ± 0.56, and 1.21 ± 0.58 mg/kg at the MW, BS, WDA, and IA sites, correspondingly, which could be generated from the metal workshop, industrial activity and Cd plated substances [59]. For Cd concentration in the soil, about 80 % and 66 % of samples exceeded TRV, and DSWS, respectively ( Table 4). The findings of this study were compared with previous studies presented in Table 5.

Geochemical indices for soil contamination assessment Contamination factor (CF), Enrichment factor (EF), Geo-accumulation index (I geo ), and Pollution load index (PLI).
The CF values for As ranged from 1.05-2.32, indicating a "moderate contamination" ( Table 6). The range of CF for Cd (2.01-6.84) displays a medium to high contamination. Pb (0.2-3.92) and Cu (0.46-4.62) were low to considerable contamination levels, whereas most of the sampling points for Cr and Ni suggested low contamination.
Conversely, it was discovered that the total CF values for all PTEs at research sites decreased the value of Cd > As > Pb > Cu > Cr > Ni, representing that Cd and As potentially create a concern for the ecosystems in the area. Considering the C d , the decreasing order of metal contamination sites as MW > IA > BF > WDA > BS > M > PS > PP > PA > PG > RA (Table 6). In the study area, PA, PG, M, PP, PS, and WDA show moderate and considerable contamination levels indicated by BF, IA, and BS. On the other hand, RA and MW show low and high degrees of contamination, respectively, due to the continuous dumping of waste materials and the discharge of industrial and urban effluent (Islam et al., 2015).
The EF is a standardized technique that is frequently used to classify the metal concentrations related to soils [15,42]. The highest value of EFs was found for Cd (24.74) and the lowest values for Cr (0.53) indicate that the studied metals were minimal to enrich soil significantly. The research locations' EFs were arranged in descending order as follows: (Table 7). For all sample sites combined, the EF of all PTEs was in the smallest to largest of Cd > As > Pb > Cu > Ni > Cr ( Table 7). The study findings indicated that anthropogenic activities (EF>1.5) are mostly responsible for the origin of PTEs in the urban soil of the Jashore district, while this study also observed lower enrichment (EF<1.5) values that show metals in the soil come from crustal origin.
For all sample sites combined, the I geo of all PTEs were in descending order as follows: Cd > As > Cu > Pb > Cr > Ni ( Table 7). The study area soils were moderate to extremely contaminated with Cd due to air emissions, battery leachates, and Cd-plated materials, while other toxic metals were uncontaminated. This outcome also matches with other studies [7,9,17,45].
PLI is the most effective technique for evaluating the degree of toxic metal pollution. PLI values in this research ranged from 0.48-2.82, with an average of 1.29 (Fig. 2). The majority of the study site soils are contaminated (PLI > 1) by PTEs, indicating that their soil quality is deteriorating gradually. Site MW (2.82) has the highest PLI values, followed by BF (1.53) and IA (1.43) (Fig. 2) sites recommended that activities from the metal workshop, industries, and emissions from the brick field make the study area more harmful for surrounding ecosystems. The results of the overall evaluation showed that the investigated area is contaminated and faces serious threats to the surrounding environment. [17] found a significant level of PLI in their study area (Jhenaidah and Kushtia district, Bangladesh).

Assessment of potential ecological risk
The results of the PER valuation revealed that the E ri value for every metal at each location was, in descending order: Cd > As > Pb > Cu Table 5 Comparison of heavy metal concentration (mg/kg) (range) in soils of Jashore districts with other studies.  (Table 8) indicating a very high potential ecological risk for their Table 6 Contamination factor, degree of contamination and contamination level of potentially toxic elements in soils collected from Jashore district, Bangladesh.  Fig. 2. PLI values in the study area.

Table 8
Potential ecological risk factor, risk index and risk level of potentially toxic elements in soils collected from Jashore district, Bangladesh.

Non-carcinogenic health risk assessment
The potential non-carcinogenic health risk of individual mental for children and adults is determined by calculating HQ. The HQ value of each metal for ingestion, dermal contact, and inhalation routes are presented in Table S6, Table S7, and Table S8, respectively. Based on USEPA [48] guideline, HQ > 1, indicates non-cancer health risk is related to over-exposure. The study result shows that non-cancer health risks associated with each metal exposure via soil ingestion, dermal contact, and inhalation were lower than one (HQ < 1), suggesting a lower health risk for children and adults. The cumulative impacts of exposed metals were computed as HI (i.e., combining soil ingestion, dermatological contact, and inhalation). The HI values for children were 0.33-1.02 with an average of 0.65 ± 0.17 as well as for adults 0.04-0.17 with an average of 0.09 ± 0.03 (Fig. 3), except for one site (MW) for children, other all sampling sites for both children and adults indicating was a lower tendency of having non-cancer health risk.
The HI value for children was 7.22 times higher than adult inhabitants, demonstrating that children may pose 7.22 times greater chances of a non-cancer health risk than an adult in the future. Moreover, the THQ value for children was 7.82, 2.74, and 2.33 times higher than adults for ingestion, dermal contact, and inhalation exposure route, respectively. The findings of the current study conform to the observations of Ajeh et al. [64], and [17]. It demonstrates that children are more vulnerable to non-cancerous health risks of studied toxic metals than adults due to their high possibility of exposing contaminated particles through touching, mouthing, and playing behaviour [64]. Mainly, the oral ingestion route and concentration of As in soil are responsible for non-carcinogenic risks to children and adults. They might suffer from diverse non-carcinogenic health hazards including vomiting, skin diseases, respiratory tract irritation, cardiovascular, intestinal pain, hair loss, and neurological problems [16]. Ingestion was the prime route of exposing toxic elements in Benin City, Nigeria also found by Ajeh et al. [64].

Carcinogenic health risk evaluation
The summary of the CRs evaluation is presented in Table 9. The decreasing order of different exposure routes for metals is ingestion > dermal contact > inhalation for both adults and children, indicating that ingestion is the most vital exposure route for creating potential cancer risk. The LTCR of studied heavy metal following descending order of As (2.10E-03) > Cd (6.69E-05) > Pb (3.25E-05) for children and As (2.74E-04) > Cd (3.17E-05) > Pb (4.33E-06) for adults ( Table 9), suggesting that As is the most concerning metal for both children and adult which may pose adverse diverge health effects.
The overall cancer risk for adults and children was found to be 3.10E-04 and 2.20E-03, respectively, which were higher than recommended values (1E-06-1E-04) ( Table 9). Table 9 also indicates that children have 7.10 times greater carcinogenic health risk than adults due to their lower body weight and higher exposure frequency. This study finding suggested that As is the main contributor to the overall cancer risks among the other studied heavy metals (Pb and Cd) in the urban soil of the Jashore district (Table 9). Though this study did not indicate any specific forms of arsenic-associated cancer in the study area. Several authors stated that long-term exposure to As might create different types of cancer such as skin, kidney, bladder, etc [16]. On the other hand, the chances of carcinogenic risk via oral exposure to metals-contaminated soil are high (Table 9), so to reduce health risks in the future, it is vital to minimize metal poisoning of urban soils in the Jashore district.

Principal component analysis (PCA), Pearson's correlation matrix (PCM), and cluster analysis (CA) for PTEs in soil samples
To understand the origin, dispersal, and mobility of the analyzed PTEs in the urban soil, multivariate investigations (PCM, CA, and PCA) are used ( Fig. 4 and Table S9). In this study, PCA, CA, and PCM were used to find studied heavy metal sources precisely. Varimax rotation was used to perform PCA based on the concentrations of the tested PTEs. About 0.709 and 157.029 (df = 15, P = 0.000) were the values for Kaiser-Meyer-Olkin and Bartlett, respectively.
The PCA results revealed those two principal elements could account for 66.092 % of the total variation in the soil and could be used to represent the heterogeneity of PTEs. The primary principal component (PC1), which contributed 48.42 % of the total variance as well as loading with As, Cd, Pb, Ni, and Cu (Fig. 4a), which displays mutual linkage in the cluster diagram (Fig. 4b) with a considerable positive correlation [ Table S9]. The contributors of As, Cd, Pb, Ni, and Cu in the study sites probably occur from road dust, the outdoor combustion of waste items, pollution from industry and automobiles, vehicle tires, excessive use of chemical fertilizers, farming practices and discharging of wastewater (e. g., rechargeable battery, metal workshops, metal plating, brake pads, arsenic sulfide, and cosmetic varnish for wood industries). The PC2 component reported 17.69 % of the overall variance. It displayed positive loading for Cr, which is separately positioned in the PCA, CA (Fig. 4a, b) and demonstrations no correlation between other metals (Table S9) signifying that the Cr had a separate lithogenic source. However multivariate analysis has some variations but it shows a good harmony for studied PTEs sources recognition. Lastly, the multivariate analysis outcome shows that manmade activities are governing rather than natural sources, which is a match with earlier findings of Bhuiyan   Fig. 3. HI value for children and adults in the study area.

Table 9
Carcinogenic risk of children and adults due to ingestion, dermal contact, and inhalation of arsenic, cadmium and lead in soils collected from Jashore district, Bangladesh. et al., [65], Islam et al., [9], and Islam et al., [5].

The positive matrix factorization (PMF) approach for metal source attribution in soil
A common factorization receptor model, PMF (version 5.0), is used to determine each toxic metal's contribution as well as to identify and evaluate the likely sources of toxic metals in urban soil. In this study, Q was reduced to control the residual matrix and the residuals that were between − 1 and + 1. This was done to verify the model and determine its applicability. This system was run 20 times, with different parameters (3 or 5) being examined each time, to achieve a decent result. The association between experiential (R 2 ) and tested value ranged from 0.82 to 0.99 indicating a strong relationship and 0.33 (Pb) is stated as the lower relationship for soil samples (Table S4). Therefore, the PMF model's distribution of the analyzed hazardous metals was accurate, and the outcomes were trustworthy. Factor 1 was vastly loaded with Pb (81.8 %) and Cu (64.7 %). Furthermore, factor 1 had a considerable amount of Ni (39.8 %), and Cd (25.7 %) (Fig. 5). Factor 1 was shown to be widely distributed throughout the research location. Lead was mainly coming from dust deposition from traffic emissions, vehicle tires, and brake pads [66]. While Cu in the study area comes from industrial activities (waste  from metal and electrical working areas), waste burning, and sewage sludge [65,67]. Hence, Factor 1 was considered a diverse anthropogenic source.
Factor 2 was dominated by Cr (78.1 %), As (40.2 %), and Cd (38.8 %) (Fig. 5). Chromium ranked 21st in the richness of elements inside the Earth's crust, causing exogenic geological processes (such as erosion and weathering) that contribute to the soils [68]. Besides, in the study area As and Cd was connected to urban and industrial pollutants, where waste from industrial activities, metal workshop, brickfield and burning coal and other fossil fuels, petrochemical might be a rational cause of As and Cd in the soils. Therefore, Cr is controlled by geogenic sources and As and Cd are anthropogenic inputs. Hence, factor 2 was attributable to diverse sources. Factor 3 was described by As (50.3 %), Ni (47.2 %), and Cd (35.5 %) (Fig. 5), where As and Cd were highly influenced by anthropogenic activities. Nickel in the soil may be resulting from the burning of sewage sludge and waste, as well as the combustion of coal, diesel, and fuel oil [69]. Finally, factor 3 is influenced by urban pollution.

Limitations of the study
This research provides important information regarding the presence of PTEs in the urban soil of Jashore, Bangladesh. In this study, seasonal variability was not taken into account. Additionally, PTEs in soil's bioaccessibility have not been determined, making it impossible to estimate actual exposure levels. Moreover, applying the risk assessment to these metrics makes the case even stronger. We defined the hazard quotient (HQ), hazard index (HI), and carcinogenic risk (CR) as the noncarcinogenic and carcinogenic risks of PTEs, respectively, where extra sources weren't taken into account. However, we examined 11 land use areas for six different toxic metals. Due to resource limitations, other metals and land-use sites could not be included in this study. These problems should be considered in future work.

Conclusions
Soil from the metal workshop, industrial area, and brick filed site were more heavily contaminated by studied PTEs than in other land-use areas, demonstrating that metal concentrations in the soils of the studied region varied with changes in land-use. Concentrations of Pb, As, Cr, Cu, Ni, and Cd in some sites exceeded background, soil quality indices, and soil quality standard values suggesting that soils were low to high polluted by PTEs. The substantial presence of toxic metals in soils was demonstrated by positive matrix factorization model fitting coefficients for PTEs, which exhibited R 2 > 0.70 and demonstrated the validity of the model's output. Mainly, anthropogenic activities were considered to enrich PTEs in the soils. According to the possible risk indices, only Cd had a very high risk for ecology, which ranged from moderate to very high-risk levels in the study area. This study indicates lower non-cancer risks of each metal for children and adults by exposing contaminated soil via ingestion, dermal contact, and inhalation. Still, these target groups had significant carcinogenic risks. Children are more suspected than adults compared to health risk vulnerabilities. However, considering the oral ingestion of PTEs, the risk to human health will increase. Therefore, proper actions have to be taken to considerably improve the soil quality and decrease the health risks to the residents. Further study is required to identify the toxico-dynamics of PTEs in the soil profile for taking decisions on human health and environmental risk management of the study area.

Ethics approval and consent to participant
This article does not contain any studies involving human participants performed by any authors. These studies evaluate the potentially toxic elements I urban soil and calculate the human health risk based on the assessed concentrations. No further toxicological experiments were conducted on any other living species. The manuscript in part or in full has not been submitted or published anywhere.

Funding
This research did not receive any specific grant from funding agencies in the public, private, or not-for-profit sectors.

Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Data Availability
Data will be made available on request.